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In thermoacoustic travelling-wave engines and other Stirling cycle devices, good performance depends 
on the material of a regenerator being in intimate contact with the gas inside it, so that each particle 
of gas oscillates in temperature following the adjacent material as it is acoustically displaced. This 
requires that the passages are small enough for temperature waves to penetrate across the gas path with 
the frequencies of interest. One type of ‘regenerator’ that is commonly used for this purpose is composed 
of multiple layers of woven stainless steel mesh, laid on top of one another in random registration. Asso¬ 
ciated with the thermal penetration is a viscous loss of pressure and this must be quantified if efficient 
engines are to be designed. 

In the literature, reliance has been placed on the correlation of steady-flow loss data for these meshes, 
but for the coarser ones operating at frequencies greater than 28 Hz, the assumption of quasi steady-flow 
is dubious and direct acoustic measurements must be made. This paper reports acoustic pressure loss 
data for meshes with 34 and 75 wires per inch taken in two configurations of impedance tube, and finds 
that the dependence on velocity is the same as in steady-flow, but that there is indeed some enhance¬ 
ment of loss for frequencies above 40 Hz. (Separation of the mesh layers is probably responsible for 
the anomalously low loss coefficients that were recorded in one set of data.) It is shown that the acoustic 
pressure losses can be correlated in terms that give the acoustic impedance more directly than the fric¬ 
tion factor correlations. 

© 2013 Published by Elsevier Ltd. 


1. Introduction 

In the performance of thermoacoustic heat engines, the charac¬ 
teristics of the material in the regenerator are of paramount impor¬ 
tance [1 ]. A temperature gradient is imposed on the regenerator by 
heat exchangers at each end and the function of the material is to 
ensure that acoustic waves passing through it are in thermal con¬ 
tact to some extent. 

In a ‘travelling wave’ engine, although the wave may be far from 
a pure travelling one, the contact should be intimate so that all the 
particles of gas follow the temperature of the adjacent material, 
and thus undergo a Stirling cycle as the gas is displaced alternately 
towards the hot end (while compressed) and the cold end (while 
expanded). This requires that the frequency is low enough for tem¬ 
perature waves to penetrate across the gas path, and hence (even 
for moderate frequencies) that the gas passages be very small. 
However, this must be achieved without incurring so much viscous 
or thermal dissipation that the work input to the acoustic wave is 
negated. Depending on the acoustic conditions, both of the mech¬ 
anisms may be important, but in a high impedance regenerator, 
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viscous dissipation tends to dominate. Thus, it is essential to know 
the acoustic pressure loss characteristics of the material so as to be 
able to quantify the viscous dissipation. 

One type of ‘regenerator’ that is commonly used in travelling 
wave devices is composed of multiple layers of woven stainless steel 
mesh, laid on top of one another in random orientation. This note de¬ 
scribes measurements of the acoustic characteristics of a 10 mm 
thickness of two different mesh layers, and compares the perfor¬ 
mance with the correlation presented by Swift and Ward [2]. The 
correlation was derived from the steady-flow data quoted by Kays 
and London [3] from the experiments of Tong and London [4], It 
has proved its relevance to thermoacoustic conditions through 
incorporation in engine calculations. However, with an acoustic 
wave there must be differential acceleration of different regions of 
the gas at the microscopic level within the mesh, so it would not 
be expected that the steady-flow data can be applied precisely. Di¬ 
rect measurements of the acoustic losses are presented here, thus 
contributing to the filling of a significant gap in the literature. 

2. Pressure loss coefficient measurements 

Schematic diagrams of two realizations of the apparatus are 
shown in Fig. 1 . They were both excited by two Pioneer sub-woofer 
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Nomenclature 



A 

face area of regenerator 

k 

thermal conductivity 

c P 

thermal capacity 

P 

acoustic pressure, or pitch of woven wire mesh 

G 

peak volumetric flow rate at regenerator mid-point 

Pa 

acoustic pressure 

K" 

dimensional loss coefficient 

n, 

hydraulic radius 

L 

regenerator depth 

t 

time, or thickness of wire in the mesh 

Leff 

effective depth of regenerator including end-corrections 

u', u 

local acoustic velocity, averaged across the passage 


for gas inertia 

Upore 

notional acoustic velocity in the pores of the mesh, 

LC 

dimensionless pressure loss coefficient 


based on porosity, at the mid-point of the regenerator 

P 

mean pressure 

X 

axial distance 

Pr 

Prandtl number 

yo 

channel half-width 

Re 

Reynolds number based on fixed length 

y. z 

transverse distances 

T 

mean temperature 

a 

thermal diffusivity 

T'aJa 

acoustic temperature, averaged across the passage 

y 

ratio of specific heat capacities 

Upore 

mean velocity in the pores of the mesh, based on poros- 


viscous, thermal penetration distances 


ity 

c 

specific acoustic impedance =z/pc 

v pk 

peak voltage driving the sub-woofer amplifiers 

1 9 

parameter in loss coefficient correlation 

z 

acoustic impedance (p/u) 

fiscale* d ve 

parameters in proposed variation of P 

C 

speed of sound 

n 

dynamic viscosity 

cficp), c 2 (<P) parameters in the correlation of Swift and Ward [2] 

V 

kinematic viscosity 

de 

equivalent diameter of the gas passage 

p 

mean density of gas 

f 

Fanning friction factor, defined by Eq. (1) 

Pa.Pa 

acoustic fluctuation in the density of the gas, averaged 

u 

Rott’s acoustic parameter obtained by integrating h v 


across the passage 


across the passage 

tdiff 

time constant of diffusion 

fk 

Rott’s acoustic parameter obtained by integrating h k 


porosity 

K(y, z) 

across the passage 

function describing the variation of the acoustic velocity 

CD 

circular frequency 


field 

Subscripts 

My. z) 

function describing the variation of the acoustic tem¬ 
perature field 

1-7 

position of microphones in Fig. 1 


loudspeakers, rated at 1 kW each, mounted in an MDF enclosure so 
that the majority of the sound emanated from a 250 mm diameter 
orifice in the side of the enclosure, in the first series of experi¬ 
ments, the sound then propagated through 530 mm of a horizontal 
300 mm diameter steel duct, halfway along which was mounted 
the mesh assembly, before exhausting to the laboratory. The 
assembly was a 10 mm sandwich of layers of plain weave mesh be¬ 
tween two supporting layers of coarse grid. These coarse grids 
were required to prevent the layers of mesh under test from sag¬ 
ging or springing apart, and they were assumed to make a negligi¬ 
ble contribution to the pressure loss. In the second series of 
experiments, suggested by one of the referees of the first draft of 
this paper, the duct was extended by 200 mm and the sandwiches 
of mesh were installed 20 mm from the open end. A second 10 mm 
sandwich of coarse mesh was placed near the orifice of the enclo¬ 
sure, to enhance the radial uniformity of the acoustic field. 


2.3. Technique for the first series 

The positions of four 6 mm holes in the 2 mm thick stainless 
steel pipe are also shown in Fig. 1 a. The tube mountings communi¬ 
cating with these holes were designed to support 13 mm B&K 
microphones. The calibration for these microphones was obtained 
before each experiment with a 1 kHz source giving 10 Pa, and as¬ 
sumed to apply at the lower frequencies of the experiment, on 
the basis of the manufacturers’ data. 

In the first series of experiments, driven by the need for high 
accuracy in the phases in particular, pressure ratios and phase rela¬ 
tionships were recorded using the same microphone. To obtain the 
phase relationships between positions, the microphone signals 
were each referenced to the current of the sub-woofer driving sig¬ 
nal. Using the ratios of magnitude and differences of phase, the 


two-microphone method [5] was then applied to deduce the 
acoustic velocities. Thus, the velocities at stations 1 and 2 for 
example were given by 


-P2+P1 cos (f 9 

ipc sin (™) 


and 


.sin (tf) 

u 2 = -1—rr^Pi 


pc 



u 1 , 


where x is the separation distance, co the circular frequency, c the 
speed of sound, and p the density of air at ambient conditions. 

Pressures and velocities at stations 1 and 2, and also 4 and 5, 
were extrapolated to the centre of the mesh, also using the pipe 
transmission equations. From these data, the pressure loss Ap 
and the velocity in the middle of the mesh, before its correction 
for the porosity to give u pore , were deduced. The velocities deduced 
by extrapolation from stations 1 and 2 were 5-20% greater than 
those from 4 and 5. The value from the loudspeaker side was taken 
to be representative in working out the loss coefficient because the 
larger phase differences there should yield higher accuracy. 

In collecting the data, the voltage generated by the amplifier 
driving the sub-woofers was fixed. The acoustic wave incident on 
the regenerator was varied by varying the frequency from 112 Hz 
down to 28 Hz, the range of interest in the present thermoacoustic 
studies. A second and higher voltage was then used to provide 
overlapping and larger values on the velocity scale. Electrical pow¬ 
ers up to 500 W were employed, but the maximum transmitted 
acoustic power in the duct was only about 1 W, because a predom¬ 
inantly standing wave was generated. Typical deduced profiles of 
acoustic pressure and velocity are shown in Fig. 2. 
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Fig. 1. (a) First configuration of the experiment with four microphone positions, (b) 
Second configuration with three active microphone positions (4, 5 and e). 


2.2. Technique for the second series 

In the second set of experiments, because of the absence of suf¬ 
ficient duct between the meshes and the open end, it was neces¬ 
sary to measure the acoustic pressure at the outlet by pressing a 
13 mm rubber bung containing a 6 mm GRAS microphone against 
the supporting grid, as shown in Fig. lb. The same design of rubber 
mounting supported two similar microphones in two of the three 
tubes at the positions shown. After checks that interchanging the 
microphones, with their slightly different phase characteristics, 



x /m 

Fig. 2. Calculated profiles of acoustic velocity and pressure along the duct from the 
open end towards the sub-woofers, showing the acceleration of the flow and the 
pressure loss through the confines of the mesh, at 80 Hz. 


did not affect the results substantially, it was decided to leave 
one microphone at position 4 and to move the second to positions 
5, 6 and the duct exit, recording the phase differences from the ref¬ 
erence microphone and normalizing the amplitude to that of the 
reference. This was necessary because of small variations in ampli¬ 
tude as the sub-woofers warmed up. Microphones 4 and 5 were 
analysed to deduce velocities and to extrapolate to the face of 
the mesh. Microphones 4 and 6 provided a check on this calcula¬ 
tion. A hot-wire anemometer probe was also installed at the exit 
to give an approximate indication of acoustic velocity there at 
the time when the amplitude of excitation was being set. 


3. Analysis of results 


The change in the real component of the acoustic pressure rel¬ 
ative to the velocity is a measure of the dissipative loss due to vis¬ 
cosity (see the Appendix A). There are also some capacitive and 
inductive effects, but these are small in comparison. In addition 
there is dissipation due to thermal effects giving rise to a change 
in the acoustic velocity. However, it is shown in the Appendix A 
that the ratio of the viscous to thermal dissipation is very large 
in these experiments, and that the change in velocity can be 
neglected. 

Anticipating the use of the data in thermoacoustic network cal¬ 
culations, an acoustic pressure loss coefficient is defined here as 
LC = real| ^^ |, the change in specific impedance, where u pore is 

the amplitude of velocity at the regenerator mid-point. This is cal¬ 
culated from the peak volumetric flow rate at that point G and the 
‘3-D porosity’ </> as u pore = 4L, where A is the face area of the mesh. 

To a good approximation, <j>= 1 — ^ sj 1 + for a mesh with 

wire diameter t and pitch p. However, the square root term which 
expresses the additional length of wire due to its tortuosity in the 
woven form is sometimes neglected. 

Of substantial interest is the relationship of the pressure losses 
in acoustic flows to those under steady-flow conditions. The stea- 
dy-flow loss data reported by Kays and London [3] have been 
coded by Swift and Ward [2] in the form of correlations of the Fan¬ 
ning friction factor with mean velocity U pore as 


/ 


2APr„ 

1KJ 


ClW>) 

Re 


+ c 2 (</>), 


( 1 ) 


where the hydraulic radius is r„ = = |-| and is V 4 of the 

equivalent diameter de. L is the regenerator thickness and the Rey¬ 
nolds number Re = '***-> The first term in this equation expresses 
the viscous component of the pressure loss, and the second, the sep¬ 
arated flow component, as illustrated by the manipulation of (1) to 
give 


A P = 


L /-lUpore 

de 2 r h 


ci(<M + 


4 L pU 2 por e 
de 2 


C2 (</>)• 


( 2 ) 


Swift and Ward [2] propose Ci ( c />) = 1268 - 3545</> + 2544</> 2 and 
c 2 (</>) = -2.82 + 10.70 - 8.6</> 2 . Defining 

y<< _ 2.fLpU pore — A P 
“ de 

the first and second terms in K" are typically of comparable magni¬ 
tude for our meshes when U pore = 3 m/s. For acoustic waves, Swift 
and Ward average the square term in Eq. (2) across the flow to 
obtain 



fiL 

K 


c i (</>) + 


pL 1 P\u pore \de 
x\ 3n fi v ' 


(4) 
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A dimensionless representation of this loss coefficient is 
K" 

LC = —, (5) 

pc 

where c is the speed of sound and, in non-isothermal situations, 
both this property and the gas density are evaluated at the mid¬ 
point of the regenerator. According to (4) and (5), LC is proportional 
to the length of the regenerator, and is a function of the Reynolds 
number and the acoustic Mach number, and of the porosity. The 
change in the specific impedance of a duct containing a transverse 
regenerator (referred to the free-stream velocity) can therefore be 
written 


Af = 


LC COLeff 

<i> j (j)C 


( 6 ) 


where the effective length L eJ y is also deduced empirically from the 
imaginary component of the pressure loss relative to the velocity. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

Pore velocity / m/s 


Fig. 3. Variation of the loss coefficient with velocity through the pores of the mesh 
for the first series. Comparison with values derived from the friction factor 
correlation of Swift and Ward [2] and the fit proposed in Section 5. The V pk values 
refer to the voltage driving the amplifiers for the sub-woofers. The frequency range 
was 28-112 Hz. 


4. Results 

The original choice of mesh (‘coarse’) for the present project 
was found in thermoacoustic experiments to have far from perfect 
thermal contact with the ~80 Hz waves of particular interest. A fi¬ 
ner one was therefore also selected, with some misgivings about 
the obviously high resistance to flow. 

The characteristic parameters of the two meshes are given in 
Table 1, together with three others on which steady flow tests have 
been conducted recently by Yu et al. [6], 

In the first series of experiments, the number of mesh layers 
was calculated on the basis of twice the wire diameter per layer 
to give a 10 mm thickness: 20 for the coarse mesh and 53 for the 
fine one. In the case of the coarse mesh, a test was also conducted 
with half the number of layers, to check the linearity of the data 
with respect to length. In the second series, 20 layers of the coarse 
mesh and 60 layers of the fine mesh were stacked. These gave a 
thickness of 9.5 mm in the centre under light compression, in both 
cases, against notional values of 10.2 mm and 11.4 mm, respec¬ 
tively, so apparently some interleaving of the wires in adjacent lay¬ 
ers is possible. The 12% smaller number of fine mesh layers in the 
same thickness of stack in the first series suggests that there was 
some springing apart. 

The results for the first series, in which a set sub-woofer voltage 
was maintained for the range of frequencies, are presented in 
Fig. 3, with a sample given in Table 1. 

The results for the second series are presented in Figs. 4 and 5 
for the coarse and fine meshes, respectively. In this case, the ampli- 
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Amplitude of Pore Velocity / m/s 


Fig. 4. Variation of the pressure loss coefficient with the amplitude of the velocity 
in the pores of the coarse mesh and the frequency of excitation in the second series. 
Comparison with the steady-flow correlation of Kays and London [3], as interpreted 
by Swift and Ward [2], 


tude of the acoustic velocity through the mesh at a fixed frequency 
was maintained approximately at a number of levels by adjusting 
the sub-woofer voltage while observing the rms hot-wire signal. 

The results are presented in terms of the pore velocity, rather 
than the Reynolds number Re = 4Up °" r)l , because Eqs. (4) and (5) 
show that the second term is actually dependent on the acoustic 
Mach number, However, this term reflects the influence of 
the ‘form drag’ of the wires in the mesh. The Reynolds number 


Table 1 

Dimensions of meshes investigated here and by Yu et al. [6], Results for loss coefficient from acoustic experiments and for steady flow, with a mean pore velocity of 3 m/s, which 
gives an acoustic displacement at 80 Hz of 6 mm. Comparison with the correlation derived by Swift and Ward [2] from the steady flow data of Tong and London [4] reported by 
Kays and London [3], 


Source 

‘Coarse’ mesh of Yu et al. 

‘Coarse’ mesh 

‘Medium’ mesh of Yu et al. 

‘Fine’ mesh 

‘Fine’ mesh of Yu et al. 

Wire diameter (mm) 

0.280 

0.254 

0.230 

0.095 

0.089 

Mesh # 

30 

34 

45 

75 

94 

Wire pitch (mm) 

0.847 

0.747 

0.564 

0.340 

0.270 

3-D porosity 

0.727 

0.733 

0.654 

0.781 

0.728 

2-D porosity 

0.448 

0.436 

0.351 

0.519 

0.450 

Hydraulic radius r h (mm) 

0.200 

0.174 

0.122 

0.084 

0.064 

Length (mm) 

50 

10 (or 5) 

50 

10 

50 

Frequency (Hz) 

20-120 

28-112 

20-120 

28-112 

20-120 

With a pore velocity of 3 m/s and 10 mm regenerator thickness at 80 Hz 





Reynolds number 

151 

132 

92 

64 

48 

Acoustic LC scaled to 10 mm 

0.01 

0.19 (1st) 

0.03 

0.44 (1st) 

0.14 



0.21 (2nd) 


0.66 (2nd) 


Yu et al. measurements giving LC in steady flow 

0.13 

- 

0.21 

- 

0.63 

Swift and Ward correlation giving LC 

0.14 

0.17 

0.32 

0.56 

0.78 
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Fig. 5. Dependence of the pressure loss coefficient with the amplitude of the 
velocity in the pores of the fine mesh and the frequency of excitation in the second 
series. Comparison with the steady-flow correlation of Kays and London [3], as 
interpreted by Swift and Ward [2]. 
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Fig. 6. Variation of the loss coefficient normalized by the steady-flow value of Kays 
and London [3] with dimensionless frequency for two meshes and two velocities. 


ranged up to 60 for the fine mesh, and up to 200 for the coarse, and 
the Mach number up to 0.009 and 0.013, respectively. What clearly 
emerges from Figs. 4 and 5 above is that there is also a dependence 
on the frequency. This can be expressed in dimensionless terms as 

dependence on 

5. Discussion 

5.1. Consistency of the results 

In the first series of results, satisfactory agreement in the scaling 
from 5 mm to 10 mm for the coarse mesh is shown, and also 
repeatability to within 2% for both meshes. However, the range 
of amplitude of the pore velocity varied with frequency. Because 
of the characteristics of the sub-woofers and enclosure, the magni¬ 
tude of the velocity rose to peaks at 40 and 80 Hz and then fell rap¬ 
idly with increasing frequency, so the associated acoustic 
displacement with the coarse mesh and 1 V pk excitation peaked 
at 40 Hz at about 8 mm, and fell very rapidly above 80 Hz, where 
it was 4 mm. Corresponding figures for 2.5 V pk excitation were 
14 mm and 7 mm, respectively, so end effects due to some of the 
interacting air particles spending significant time outside the mesh 
might have affected the data for this case, but the agreement of the 
data for 5 mm and 10 mm lengths suggests otherwise. 

There remained the possibility that an effect of frequency was 
superimposed upon that of velocity. However, since the results 
for different intensities of excitation (1 and 2.5 V pk in the sub-woo¬ 
fer driving signal) agreed adequately at the same velocity, but dif¬ 
ferent frequency, this possibility seemed to be negated. There was 
a wide enough range of data to establish a closely linear depen¬ 
dence of the loss coefficient on the velocity that was independent 
of frequency. 

However, when the second series of experiments was under¬ 
taken, one objective was to look systematically for any influence 
of frequency. This emerged in the data for both the coarse and 
the fine meshes. The results for all frequencies, apart from 
100 Hz with the fine mesh (Figs. 3 and 4), agreed closely with 
the Kays and London correlation [3] for steady flow at low acoustic 
velocities. However, as the amplitude increased, the higher fre¬ 
quency loss coefficients became greater than the steady-flow ones. 
Unfortunately, in the case of the coarse mesh, the results for 40 Hz 
fell below the steady correlation, so that when the enhancement in 
loss coefficient over the correlated values is plotted against the 
dimensionless frequency KM (Fig. 6), the two meshes show sim¬ 
ilar variations individually but an inconsistent variation when 


combined. Because the penetration depth is becoming comparable 
to the hydraulic radius, a greater enhancement would be expected 
in the coarse mesh case. An increase in the coarse mesh loss coef¬ 
ficients by about 30% relative to the steady-flow values would be 
required to fully rationalize the data. 

The other objective of the second series was to check that the 
thermal dissipation, which was reduced at the new location of 
the mesh in the duct, was indeed negligible. Substitution of the 
measured impedances on the mesh faces into the expression 
(A20) in the Appendix A suggested that the changes in acoustic 
velocity were insignificant compared with those in the pressure 
in the first series, and even smaller in the second series. 

The correspondence of the first and second series of data for the 
coarse mesh (two assemblies, possibly of slightly different thick¬ 
ness) was to within the estimated precision of measurement (see 
Section 5.2), given the range of frequencies superimposed in the 
first set. However, there is a clear discrepancy between first and 
second sets in the fine mesh results, the nominally 12% thicker sec¬ 
ond stack giving losses some 30% greater than the first and agree¬ 
ing closely with the steady-flow correlation at low frequency. 
Other tests have indicated that this is most likely explained by 
the poor compaction of the stack in the first case. 

5.2. Assessment of accuracy 

Clearly, adjacent pairs of holes are closer (0.33 pipe diameters) 
than is desirable for the two-microphone method because of the 
requirement for very accurate phase measurement. However, the 
shortest wavelengths were at least ten times this diameter so the 
wave should have been planar to a high degree. As a check on 
the accuracy of the data, a sensitivity analysis was carried out to 
assess the possible effect of systematic errors in the pressures on 
LC, particularly in recognition of this very short distance between 
the pressure measurement locations. 

Comparison of the magnitude of the velocities deduced for po¬ 
sition 4 in the first series from the pairs of microphones on the two 
sides of the mesh showed discrepancies ranging from 1% to 12%. 
Conceivable errors in the pressure magnitudes could not account 
for the discrepancies in velocities, which were predominantly in 
the real components. However, they could have been accounted 
for by an error of -0.8° in the phase of either p 2 or p 4 . A systematic 
measurement error of this magnitude is possible, for instance 
because of the positioning of the microphones relative to acoustic 
waves that were not truly planar at all locations. 

Systematic errors in the pressure magnitudes at particular loca¬ 
tions would have had a significant effect on the deduced values of 
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loss coefficient. Typically, a 1% decrease in Ip^ or a 1% increase in 
|p 2 | gave rise to a 7% increase in LC, while changes in |p 4 | had much 
less influence. The possible error of -0.8° in the phase of p 4 de¬ 
duced above would have given rise to the measured LC being 
10-20% too low. Repeatability suggests random errors of up to 
0.4° in the phases, and the resulting 10% error in LC is consistent 
with the scatter in Fig. 3. 

In the second series for the coarse mesh, velocities and pres¬ 
sures at the mesh location extrapolated from the microphone pairs, 
|p 4 | and |p 5 |, and |p 4 | and |p 6 |, were within 10% of one another at 
40 Hz and 60 Hz, but often differed by more at the higher frequen¬ 
cies. However, the boundary layer perturbation at the join between 
the two duct sections may have affected the accuracy of the data 
based on |p 4 | and |p 6 |, so they were not used in the final analysis. 

These are clearly only indications of the sensitivity of the results 
for LC, and the systematic uncertainties in the measured parame¬ 
ters cannot be precisely established. Overall, systematic errors of 
±10% must be allowed and this goes some way to explaining the 
inconsistency in the coarse and fine mesh results in Fig. 6. 

5.3. Comparison with other data 

Wakeland and Keolian [7] investigated the friction losses of a 
wide variety of single-mesh screens and showed that / was well 
correlated by a function of the ‘orthogonal’ or ‘2-D projected’ 
porosity and by a Reynolds number based on r h . However, their 
maximum frequency was only 5 Hz, so on di ff <c 1, where T di ^_ ", 
the characteristic time for viscous diffusion. Since this is essentially 
the same criterion as ^ <r 1, quasi-steady behaviour would be ex¬ 
pected. This is also true of the oscillating flow data for stacked 
screens obtained by Tanaka et al. [8] (up to 10 Hz) and by Hsu 
et al. [9] (up to 4 Hz). The latter did indeed deduce that steady flow 
correlations applied. 

At higher frequencies, some recent data for steady and acoustic 
flow through meshes are given by Yu et al. [6], Results read directly 
from their graphs (with very limited extrapolation and not relying 
on their correlations) are compared with the present acoustic re¬ 
sults in Table 1 for a pore velocity of 3 m/s. Their data are for 
30#, 45#, 94#, 180# and 200# meshes. The first three are desig¬ 
nated ‘coarse’, ‘medium’ and ‘fine’ to correspond approximately 
with the two meshes investigated in the present experiments. Re¬ 
sults for the two even finer meshes are not included in the table. 

The steady flow friction factor data of Yu et al. [6] are somewhat 
lower than the Swift and Ward 2] correlation if values for porosity 
are calculated from the published wire dimensions using the 
expression in Section 3. This is particularly true of their ‘medium’ 
and ‘fine’ meshes, which are lower by 35% and 20%, respectively. 
However, the acoustic data in [6] are lower by almost an order of 
magnitude, so these data are certainly not consistent with the 
present results. 

5.4. Alternative formulation of viscous loss 

From the analysis in the Appendix A, it may be postulated that a 
reasonable second-order representation of f v for the tortuous paths 
represented by randomly registered meshes is 

(7) 

where 0 is an empirical function of the Reynolds number, and pos¬ 
sibly the non-dimensional frequency, ^ for each type of 

mesh. It is seen that 


which is indeed independent of frequency if 9 is. Moreover, if 0 has 
the functional form 


dscale 


1 T | Upore 


(9) 


then LC is linearly dependent on velocity, as observed. 

Taking this prescription, the best fits shown in Fig. 3 are 

r 1 = 1.25 + 0.0048 Re, (10) 

for the coarse mesh, 

r 1 =2.08+ 0.0237 Re, (11) 


for the fine mesh. 

The original formulations of (13) and (14) were based on the 
assumption that the velocity through the holes of the mesh would 
be an important correlating parameter, and hence the losses 
should depend on the 2-D projected porosity. Whichever porosity 
is selected, the coefficients in these correlations are very different 
for the two meshes, so no universal prescription can be proposed. 

As seen in Fig. 3, very satisfactory, correlations have been estab¬ 
lished to describe the viscous characteristics of the two meshes in 
the first series of tests over a wide frequency range in a form that is 
useful in thermoacoustic calculations, to a precision of ±10%. How¬ 
ever, the results in the second series of tests indicate that much of 
this scatter can be eliminated by acknowledging that 6 is actually 
frequency-dependent. Moreover, the difference in the fine mesh 
results from the two tests suggests that the tightness of packing 
is also a factor. 


6. Conclusion 

Acoustic pressure losses have been measured for two grades of 
randomly-orientated stainless steel, woven mesh (34# and 75#) in 
two different configurations over the frequency range 28-112 Hz. 
The change in impedance in all four cases has been shown to have 
a component that is linearly-related to the acoustic velocity. Three 
of the data sets converge at low amplitude and low frequency 
within the expected accuracy of the measurements to the values 
derived by Swift and Ward [2] from the steady-flow correlations 
of Kays and London [3 ], based on the data of Tong and London 
[4], However, the fourth set (for the 75# mesh) have 20-30% lower 
pressure loss than predicted by Swift and Ward [2], possibly due to 
poor compaction of the meshes. 

For the results from the second configuration, a correlation of 
the departures from the steady-flow values with frequency (made 
dimensionless by the viscous penetration time) was attempted. 
The individual data sets both increased monotonically with fre¬ 
quency. Geometrical differences may have been responsible for 
causing these coarse mesh data to lie some 25% below where they 
were expected from the fine mesh results. 

Taking the formulation of the momentum equation developed 
extensively by Swift [ 1 ], the independence of frequency in the loss 
data has been shown in Appendix A to be consistent with the solu¬ 
tions for a variety of straight-channel types of regenerator in the 
low-frequency limit of large viscous penetration depth. This leads 
to a simple empirical correlation of the data for pressure loss, with¬ 
out invoking friction factor correlations, although they are equiva¬ 
lent. It is also shown that the acoustic pressure losses are much 
greater than the acoustic velocity reductions due to thermal dissi¬ 
pation, on a percentage basis. 
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Appendix A. Viscous and thermal losses in acoustics 

In Rott’s acoustic approximation [10], as developed by Swift [1], 
the momentum equation can be written 


du' dp„ 


P 


d 2 u' tfu 


| 2 ,,/ 


dy 2 dz 2 

with u' = 0 at a solid surface. This can be transformed into 
du' m u , dp a 

-My,*)]-£-, 

and then the unknown function h v (y, z ) satisfying 

d 2 h v d 2 h v _ n .h v 

~W + ~d?~ ft 

integrated with respect to y and z to obtain 
dp a _ imp 
dx (1 — fv ) U ’ 


(Al) 


(A2) 


(A3) 


(A4) 


where p a and u are the one-dimensional parameters, averaged over 
the transverse directions, and f v is obtained by averaging h v in those 
directions. Thus the viscous pressure loss across an obstruction can 
be written 


A p a icoAx 


(A5) 


Similarly the energy equation in the absence of an axial temper¬ 
ature gradient, 


pC Oh — _)_ k [ if 'a _(_ “ 

P p dt dc + \ dy 2 a * 2 


dz 2 


can be transformed into 

pCpT' a = [1 - h k (y,z)]p a , 
where 


(A6) 


(A7) 


Swift [ 1 ] shows that the functions/are very similar for a variety 
of straight channels when r h IS Vik < 0.8, where 5 V = and 

S k = yS are the viscous and thermal penetration depths, respec¬ 
tively. For parallel plates, the accurate expression is 
tanh{(l +i)f} 


/ = 


(A14) 
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where y 0 is the channel half-width. When expanded, it can be 
shown that in the limit r h /<5 = y 0 /<5 -c 1 
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We may postulate that a reasonable second-order representa¬ 
tion of f for the tortuous paths represented by randomly-oriented 
meshes is 


'-’-‘ Is . 


8, 


(A16) 


where 0 is an empirical function of the Reynolds number for each 
type of mesh. Thus, defining the non-dimensional pressure loss in 
Eq. (A5) as LC, 

2 ‘ r ii /i i; 

(A17) 
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which could be independent of frequency if 9 is. 

Moreover, since f k is obtained by an integration of the solution 
to the same equation with the same boundary condition as/ v , apart 
from the change in 5, we may take 

(A18) 


* = 1 -'£)*■ 


with the same parameter 0. The ratio of the two contributions to the 
dissipation in Eq. (A13) is 

viscous dissipation (pc|u[ 


thermal dissipation 
Moreover, 


|Pa I 
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d 2 h k 0 2 h k _ n .h k 
dy 2 + dz 2 sl ’ 


(A8) 


with T' a = 0 at an isothermal surface. Then, invoking continuity, 


imp a +-^(P u ) = 0. 

and the perfect gas relationship 

Pa = Pa,T 1 

P p~ T ' 


(A9) 


(A10) 


du 


(All) 


so that the velocity change due to thermal relaxation can be written 
icoAx, 


pcAu = - [1 + (y - 1 )f k ]p a . 


(A12) 


Thus the total acoustic dissipation associated with the obstruc¬ 
tion is 


-^Re[A(p a u)] = ^Re[p a Au + Ap 0 u] 


Aoo Ax 

^c~ 


Im 


{1 + (y - l)/ fc }if^2_ + 


I Pal , PCW 
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(A13) 


and for the values measured in the first series, this ratio is <0.02 at 
100 Hz and <0.01 at 40 Hz. Thus the pressure loss is far more signif¬ 
icant than the change in velocity in these circumstances. 
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